XXZ spin-1/2 representation of a finite-(7 Bose-Hubbard chain at half-integer filling 
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Using a similarity Hamiltonian renormalization procedure, we determine an effective spin-1/2 
representation of the Bose-Hubbard model at half-integer filling and at a finite on-site interaction 
energy U. By means of bosonization, we are able to recast the effective Hamiltonian as that of a 
spin-1/2 XXZ magnetic chain with pertinently renormalized coupling and anisotropy parameters. 
We use this mapping to provide analytical estimates of the correlation functions of the Bose-Hubbard 
model. We then compare such results with those based on DMRG numerical simulations of the Bose- 
Hubbard model for various values of U and for a number L of lattice sites as low as L ~ 30. We find 
an excellent agreement up to 10% between the output of analytical and numerical computations, 
even for relatively small values of U. Our analysis implies that, also at finite U, the ID Bose-Hubbard 
model with suitably chosen parameters may be seen as a quantum simulator of the XXZ chain. 
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I. INTRODUCTION 

The study of magnetic systems is one of the more active 
fields of research in condensed matter physics 1 : the vari- 
ety of emerging ground-states, as well as the rich phase 
diagram of magnetic lattices, makes these systems an 
optimal testbed where it is possible to study the com- 
petition between various orders and frustration effects 2 . 
From this perspective, it would be very useful to be able 
to engineer synthetic physical systems effectively describ- 
ing magnetic model Hamiltonians, with tunable geometry 
and parameters. 

A promising route is provided by cold atoms setups: 
for instance, itinerant magnetism in bulk ultracold Fermi 
systems with repulsive interactions has been experimen- 
tally studied 3 , while small spin networks have been simu- 
lated with ion chains 4 . Effective nearest-neighbour spin- 
spin interactions for atoms in neighbour wells of an op- 
tical lattice may result from super-exchange couplings: 
the corresponding second-order tunneling has been ob- 
served in array of double wells 5 . Furthermore, using fast 
oscillations of the optical lattice, one can control the sign 
of the nearest-neighbour tunneling 6 , which has been re- 
cently used to simulate classical frustrated magnetism 
in triangular lattices 7 . Finally, in ultracold atomic sys- 
tems one may use two-component gases: here, the two 
internal degrees of freedom correspond to the simulated 
(pseudo)spins. In optical lattices, in the regime of strong 
interactions, by adjusting the external potential one can 
control spin interactions between spin states of neigh- 
bouring atoms 8 . The recent experimental realization of 
controllable Bose-Bose mixtures 9 then paves the way to- 
wards the experimental simulation of spin Hamiltonians, 
in which the atomic counterpart of magnetic phases, like 
the antiferromagnctic Neel and the XY ferromagnetic 



phases (corresponding respectively to the checkerboard 
and the supcrcountcrfluid phases 10 ) may be detected and 
studied. 

A key tool in the manipulation of ultracold atomic sys- 
tems is the possibility to superimpose and control optical 
lattices 11 . The low-energy properties of ultracold bosons 
in deep optical lattices are well captured by the Bose- 
Hubbard (BH) Hamiltonian 12 
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In Eq. (1), i, j denote the lattice sites, while stands 
for any pair of nearest neighbouring sites. The operators 
b\ (bi), with [6^ , &t] = Sij and n.; = b\bi, create (annihi- 
late) a boson in the site i. The parameter t denotes the 
hopping strength, and U (V) is the interaction energy of 
two particles at the same site (at two nearest neighbour- 
ing sites). 

The use of optical lattices in ultracold atomic sys- 
tems is also central in other proposals to simulate spin 
Hamiltonians: in Rcf. 13 it was suggested to use spinor 
Bose gases in optical lattices to implement the quadratic- 
biquadratic spin Hamiltonian. The mapping of a Mott 
insulator of one-component ultracold bosonic gases in a 
tilted ID optical lattice onto a spin Hamiltonian was dis- 
cussed in Ref. 14. Based on this proposal, the quantum 
simulation of an antiferromagnctic spin chain (more pre- 
cisely, an Ising model in a transverse field 15 ) has been ex- 
perimentally demonstrated 16 : varying the effective mag- 
netic field, a quantum phase transition occurs between a 
paramagnetic phase and an antiferromagnetic one. These 
phases (and the corresponding phase transition) were de- 
tected by measuring the probability to have an odd occu- 
pation of sites, while the formation of magnetic domains 
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was observed using in-situ site-resolved imaging and noise 
correlation measurements 16 . 

Apart from the simulation of quantum magnetism, 
more generally the possibility of realizing BH models 17 
through the use of ultracold atoms in optical lattices 12 
attracted a huge interest in the last decade for a num- 
ber of other reasons: the experimental detection of the 
Mott-supcrfluid transition 18 motivated a thorough study 
of the quantum phases of the BH model. From this point 
of view, arrays of ultracold atoms are versatile synthetic 
systems for the study of strongly interacting lattice sys- 
tems 19 ' 20 in which it is possible to investigate a huge va- 
riety of phenomena, ranging from quantum phase tran- 
sitions 11 to superfluid dynamics in lattices 21 and non- 
equilibrium dynamics of quantum systems 22 . 

For very large values of U, i.e. for t/U <C 1, the BH 
model may be mapped into the Heisenbcrg XXZ spin-1/2 
Hamiltonian 

#xxz = J Y, + s i s ) - As ? s l) . ( 2 ) 

(i,3) 

where s,; = Si/ 2 = (sf, s?, sf) are the 5 = 1/2 spin 
operators, Si being the Pauli matrices, J is the nearest- 
neighbour coupling, and A is the anisotropy parameter 
(A = ±1 respectively correspond to the antiferromag- 
netic and the ferromagnetic isotropic Heisenbcrg model). 

The use of lattice spin systems for interacting bosons 
traces back to the classical papers by Matsubara and 
Matsuda in the 50's 23 . In those papers the properties of 
helium II are studied assuming that each atom in liquid 
helium can occupy one of the lattice points: this assump- 
tion naturally leads to a (BH-like) lattice Hamiltonian. 
The further assumption that two atoms cannot simulta- 
neously occupy the same lattice site (due to the hard- 
core part of the interparticle interaction between Helium 
atoms 24 or, in BH language, by a an infinite U) leads to 
a XXZ spin-1/2 model in a magnetic field 23 . To quali- 
tatively understand the emergence of a spin representa- 
tion of the one-component BH model one may say that, 
when U — > oo and if two states per site give a dominant 
contribution to the energy, an XXZ Hamiltonian is re- 
trieved: this is exactly what happens when the filling /, 
defined as the average number of bosons per lattice site, 
is half-integer. Indeed, when f = n + 1/2, with n inte- 
ger, the relevant states in the Fock space are given by |n) 
and \n + 1) (deviations from half-integer fillings would 
result in a magnetic term in the XXZ Hamiltonian). For 
half-integer /, at the leading order in t/U —> one has 
J = 2t(f + 1/2) and A = V/J (see the discussion in 
Section III). 

The XXZ model is a paradigmatic spin Hamilto- 
nian which has been the object of many investiga- 
tions (see e.g. Ref. 25) and that in ID is exactly solv- 
able by Bethe ansatz 26-28 ; this provides an ideal arena 
to test different analytical and numerical techniques, 
from bosonization 29,30 to density matrix renormalization 
group (DMRG) 31 . The study of (static and dynamical) 
correlation functions in this model is currently an active 



area of research 32-41 and exact analytical results for the 
correlation functions at small distance (both at zero and 
finite temperature) are by now available 34 . The asymp- 
totic form of the ground-state correlation functions in 
the thermodynamic limit is power-law with an exponent 
that has been obtained by comparing the result of abelian 
bosonization with the Bethe ansatz solution 42 : for an 
open chain in the region — 1 < A < 1, the numerical find- 
ings for correlation functions obtained with DMRG were 
compared with the results of a low-energy field theory, 
showing a very good agreement and allowing for precise 
estimates of the amplitudes of the correlation functions 43 . 
In turn, the obtained amplitudes were found in agree- 
ment with the analytical expressions given by Lukyanov 
and Zamolodchikov 44,45 . Finally, exact results obtained 
for the XXZ chain in a special scaling limit were used 
to compute the local correlations of a continuous Lieb- 
Liniger ID Bose-gas 46 . 

In this paper we determine a correspondence between 
the BH chain at half-integer filling for finite U and a ID 
XXZ spin-1/2 model. This enables us to provide analyt- 
ical expressions for the BH correlation functions, which 
we compare with numerical results obtained with DMRG, 
showing that there is a very good agreement both at large 
and small distances and also for U/J as low as ~ 2 and 
for a number of sites L > 30. As a consequence, the nu- 
merical determination of the phase transitions superfluid- 
charge density wave and supcrfluid-Mott insulator (re- 
spectively corresponding, in the effective XXZ chain, to 
A c ff = 1 and A e g = — 1) well agrees with the analytical 
results for the XXZ chain. Besides providing a way to an 
analytical study of the quantum phase transitions of the 
BH model, our analysis may reveal very useful in deter- 
mining the structure factor as well as the response of the 
BH chain to external perturbations. As a result, we are 
able not only to provide analytical expressions for ID BH 
correlation functions, but also to show that by realizing 
the ID BH chain at half-integer filling one can provide a 
quantum simulation of the XXZ chain. 

In the following we shall derive an effective spin-1/2 
Hamiltonian for the BH chain at half-integer filling as 
a power series of t/U. Following Refs. 47,48, we per- 
form a continuous unitary transformation S which block- 
diagonalizes the BH Hamiltonian in the basis of the eigen- 
vectors of _ffgfj with t = 0. The latter requirements yields 
an equation for S which we solve perturbatively to the 
order (t/U) 2 . We finally show that, using bosonization, 
this Hamiltonian can be recastcd in the XXZ form with 
pertinent coupling and anisotropy parameters. We ob- 
serve that, while to the first order in t/U one finds a 
XXZ model with J = 2t(f + 1/2) and A = V/J, to the 
next order in t/U one gets an effective spin Hamiltonian 
which is not of the XXZ form, since it contains also next- 
nearest neighbours and 3-spin terms (this is the bosonic 
counterpart of a similar computation done for the 2D 
Fermi-Hubbard model 49,50 , where 4-spin terms appear). 
However, in ID it is possible to further proceed using 
bosonization: we do this for the BH chain, where one 
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can obtain analytical results introducing a Luttinger liq- 
uid description of the effective Hamiltonian enabling to 
incorporate the long-wavelength behaviour of non-XXZ 
terms in the effective coupling and anisotropy parame- 
ters, J e fi and A e ff, which are now function of t, V, f and 
U. 

The plan of the paper is the following: in Section II we 
introduce the BH and the XXZ models, recalling some 
useful properties and results. In Section III we employ the 
continuous unitary transformation introduced by Glazek 
and Wilson 47 to approximate the BH chain at half-integer 
filling with an effective spin- 1/2 effective Hamiltonian. 
In Section IV we use bosonization to recast this effective 
Hamiltonian as an XXZ Hamiltonian, with effective cou- 
pling J e ff and anisotropy A c ff. In Section V we establish 
the correspondence between the correlation functions of 
the BH model and the ones of the XXZ chain. The results 
of Section III-V are then used to compare the analytical 
results obtained for the BH chain correlation functions 
with the numerical findings obtained by DMRG numeri- 
cal simulations: the results of this comparison are detailed 
in Section VI, both for the correlation functions and the 
phase transition points. Section VII is devoted to our 
conclusions, while more technical details are contained 
in the Appendices. 

II. MODEL HAMILTONIANS 

In this Section we review the basic properties of the 
BH and of the spin-1/2 XXZ Hamiltonians, in particular 
focusing on known analytical results about the real space 
spin correlations in the XXZ chain. 

A. Bose-Hubbard model 

The low-energy properties of interacting bosons in a 
one-dimensional deep optical lattice are in general well 
described by the Bose-Hubbard Hamiltonian (1), which, 
in ID and with open boundaries, reads: 

L-1 L 

i=l i=l 
L-1 

+ V n » n *+i- ( 3 ) 

i=l 

We denote with N the total number of particles in the 
L-site chain, so that the filling /, that is, the average 
number of particles per site, is given by / = -j- . For al- 
kali atoms usually V <C U, but with dipolar gases (or 
polar molecules) V could be comparable with U: exper- 
iments with dipolar gases 51 and long-lived ground-state 
polar molecules 52 in optical lattices have been already 
performed (see also the review 53). 

A large amount of experiments investigated the prop- 
erties of the BH model: the main reason for this interest 



lies on the fact that this model exhibits a quantum phase 
transition between a superfluid phase (for t/U 3> 1) and 
a Mott insulator (for t/U < l) 17 . A finite V gener- 
ally favours charge density wave phases: e.g., for half- 
integer filling / = 1/2, a large V S> t, U will result in a 
ground-state of the type |1, 0, 1, 0, • • ■ ) (where in general 
|ni, ri2, ri3, ■ ■ ■ ) is an eigenfunction of i?BH with t = 0). 
The ground-state of the BH model has been studied in 
the seminal paper in Rcf. 17 using the grand-canonical 
ensemble, where the chemical potential [i is introduced 
to enforce the constraint on the number of particles. The 
phase diagram in the U — \i plane shows the characteristic 
lobes 17 : for a pertinently fixed value of fi the half-integer 
fillings may be made correspond to the "basis" of the 
lobes (i.e. where the lobes touch) and, for V = 0, one has 
a superfluid for each finite value of t, while a finite and 
positive value of V gives rise to a charge density wave 
region among the Mott lobes. 

The Mott-insulator/supcrfluid transition was first ob- 
served in 3D 18 and subsequently in ID 54 and 2D 55 . The 
effect of a superimposed external potential (typically a 
parabolic one) has been also considered: the so-called 
wedding-cake-like density has been studied both theo- 
retically 56,57 and experimentally 58,59 . The phase coher- 
ence properties of ultracold bosons in optical lattices have 
been studied, as well, showing that phase coherence on 
short length scales still persists deep in the insulating 
phase 60 . The BH model in a ID geometry can be ob- 
tained either by tightly confining the bosonic cloud in 
two radial directions in presence of a periodic potential 
in the transverse direction, or by creating many (even- 
tually uncoupled) tubes with a 2D optical lattice. The 
Bragg spectroscopy of interacting one-dimensional Bosc 
gases loaded in an optical lattice across the superfluid 
to Mott-insulator phase transition has been performed 
and the properties of the strongly correlated phases in- 
vestigated looking at the Bragg spectra 61 . The excita- 
tion spectrum in the strongly interacting regime has been 
studied also in presence of a tunable disorder, created by 
a bichromatic optical lattice, showing a broadening of the 
Mott-insulator resonances 62 . 

The finite-V ID BH model has been studied with a 
number of analytical and numerical techniques: in par- 
ticular in Rcf. 63 the phase boundaries of the Mott insu- 
lators and charge density wave phases were determined 
by DMRG. The zero-temperature phase diagram both of 
the BH model and of a spin-S Heisenbcrg model was con- 
structed and their relation investigated 64 . The role of V 
in inducing supersolid phases in the ID BH model was 
also studied 65-68 . Bosonization techniques quite generally 
provide in ID a very effective way to compute the corre- 
lation functions and their decay at large distance: thus, 
they have been applied as well to ID BH chains (see the 
review 69). 

Finally, we mention that the effect of intersite inter- 
actions was considered since the 90's in a related model, 
the quantum phase model, describing Josephson junc- 
tion arrays 70 : the quantum phase model can be obtained 
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from the BH model for large filling per site when the 
number fluctuations are negligible in the kinetic term. In 
the correspondence between the two models the chemi- 
cal potential \i term in the BH model corresponds to the 
so-called "offset charge" q, which are external charges 
present in the superconducting network 70 : the lobes in 
the quantum phase model are equal, since there is an 
invariancc for q — > q + 2e (2e being the charge of the 
Cooper pairs), and an half- integer value of the filling / 
corresponds to half-integer values of the offset charges 
q/2e. The study of intersite interactions is relevant in 
Josephson junction arrays since the interaction term de- 
pends on the capacitance matrix C^, which is in gen- 
eral not diagonal, resulting in terms of the form VijUiUj, 
where Vij oc (Cy) : as a mean-field analysis shows 71 , 
for a diagonal capacitance matrix one has that at T = 
the superconducting phase is obtained for each value of 
the Josephson energy Ej (oc t in the mapping) and that 
at q = e one has a finite critical temperature for the 
Mott-insulator/superfluid transition for each finite value 
of Ej (unlike q = 0, where a critical value of Ej is re- 
quired). Non-diagonal terms of the capacitance matrix 
favour charge density waves 70 : the role of the intersite 
terms was considered for superconducting chains and the 
corresponding phase diagram investigated 72,73 , revealing 
that in ID a (superconducting) repulsive Luttinger liquid 
phase exists. 

To conclude this Section let us mention that, in the 
rest of the paper, we will mostly deal with half-integer 
fillings, / = n + i, with n = 0, f , 2, • • • . The reason for 
such a choice is that in this case the relevant states for 
the description of system for U — > oo are just |n) and 
|n + f). Simple arguments - reviewed in Section II B - 
then show that, to first order in t/U, the BH Hamilto- 
nian is mapped into an XXZ spin-I/2 Hamiltonian which 
is integrable in ID. Within the XXZ-modcl framework, 
it is also possible to consider small deviations from the 
half-filled regime, which mainly give rise to a uniform 
magnetic field in the z-direction. Even though we will not 
consider large fluctuations in / (of order 1), it is possible 
to take them into account, by keeping, as relevant states 
for U — ^ oo, |n), \n— 1), \n+ 1). In this case, an effective 
spin-1 XXZ effective model (in general not integrable) 
is expected 74 . Spin- 1 models exhibit a gapped (Haldanc) 
insulator phase 75,76 , which has been investigated in the 
context of the ID BH model 77 " 81 . 



B. XXZ chain 

For a chain with L sites and open boundaries, the 
Hamiltonian of a spin-1/2 XXZ model given in Eq. (2) 
particularizes to: 

L-l 

H XXZ = ~JJ2 ( s ^+i + s * s *+i - As ^+i) • ( 4 ) 
i=l 



The global minus sign in the couplings has been intro- 
duced in order to more easily perform the comparison 
with the BH model, and it can be readily gauged away 
by implementing the canonical mapping to the spin-1/2 
operators r" defined as t//' v = (— l) J s^' y , t* = s|. There- 
fore the chain is antiferromagnetic (ferromagnetic) for A 
positive (negative). 

Following Ref. 23, one can derive the Hamiltonian in 
Eq. (4) from the BH Hamiltonian (3) at half-integer fill- 
ing / and for U — >■ oo. To do so, let us define = Uj — / 
(so that the eigenvalues of s* are ±5). Since for t = the 
energy per particle is (for L — > 00) e = Uf(f—l)/2+Vf 2 , 
it follows that i?BH — > VsjSj +1 , i.e. 

J A = V. (5) 

Similarly, for / ^> 1, one gets J ss 2tf as one can see by 
putting bi ~ \f]e l '^ i and mapping the obtained result in 
the XXZ spin-1/2 language 72 : for finite values of / one 
gets (see Section III) 

J = 2t (/ + 5) ' (6) 

Eqs. (5) and (6) provide the desired mapping between 
the BH model and the XXZ Hamiltonian to lowest order 
in t/U. However, as we are going to see in Section V, 
to get a quantitative agreement between the BH and the 
XXZ correlation functions even for t/U relatively small 
(as low as 0.1 for / = 1/2) one has to go to the next order 
in t/U: the corresponding Hamiltonian is determined in 
Section III and recasted in XXZ form via a Luttinger 
representation in Section IV. We remark that, since our 
result are obtained at half-integer filling, we may omit the 
addition of a magnetic field term of the form oc s| 
to Eq. (4). Indeed such a term is proportional to the 
total spin = s i ln ^ ne z direction and, since the 

system is half-filled, only eigenstates of Hxxz with Sj, = 
are physically meaningful - notice that in the following 
analytical results based on the XXZ Hamiltonian (4) are 
compared with numerical DMRG simulations of the BH 
chain in the canonical ensemble, where ^ ■ rij is conserved 
and equal to N. 

The Hamiltonian i^xxz is exactly solvable by means of 
standard Bethe ansatz techniques 26-28 : however, explic- 
itly computing the real space spin-spin correlation func- 
tions is quite a difficult task. Exact analytical results for 
short-range correlators in a range of up to seven lattice 
sites were reported for the isotropic Heisenberg model 
in Ref. 40, in the thermodynamic limit (L — > 00) and 
at arbitrary finite temperature, and for finite chains of 
arbitrary length L in the ground-state. Results for short- 
range correlation functions are also available for the XXZ 
chain 34 . For large distances, using the standard bosoniza- 
tion approach 29,30 to spin-1/2 XXZ model 82 , one may 
find out all the spin-spin correlation functions in terms of 
two-point correlators of pertinent conformal operators 43 : 
in the thermodynamic limit one finds the asymptotic 



forms 



A, 



^ois?4i^o) = (-irv 



-4, 



A, 



i/n ' (7) 



|*-i|" |i-j|" +1 / r ' 
where |^o) is the ground-state of flxxz and we set 



7] = 1 arccos A. 

7T 



(8) 



Analytical expressions for the correlation amplitudes 
A x , A x and A z entering Eqs. (7)-(7) were presented in 
Refs. 44,45 and further discussed in Ref. 83 (see also the 
discussion in Section V of Ref. 41): 



A x 
A, 
A T 



with 



dt 



8(1-77)2 

2/7(1 - V) 
2A l ' r > T 



sinh^t) 



(9) 
(10) 

(11) 



t \ sinh(i) cosh[(l — rj)t] 



cosh(2?7i)e~ 



1 



t \ 2sinh(?^) sinh(i) cosh[(l — r])t] 



+ 



sinh(?/i) r\ 



V 2 + i e _* 



dt ( sinh[(2?7 - l)t] 2r] - 1 

t \ sinh(?7i) cosh[(l — rj)t] r] 



-21 



and 



A = 



(2(1%)) 



(12) 



and r(x) being the Euler's Gamma function. 

Analytical expressions (in the large-L limit) for the 
subsequent prefactors of the correlation functions are re- 
ported in Refs. 36,41. 

For chains of finite size L with open boundary condi- 
tions, one obtains 43 : 



and 



(M s i s i IV'o) 



2/i (2i)f 1 (2j) \h(i- j) h (i + 5) / 4^77 \f 2 (i- j) / a (i + j) 



a ( (~l) i 

27TT7 \ / 1 (2i) 



b(* - J) + 9{i + j)} - 



(-1) J 

/^(2i) 

2?7 



(V>o|s?4l^o) = 



/ f (2»)/ f (2j) 
fr,(i~ j)f v {i + j) 



[fif(« - i) - + i)] 



be 



sgn(i - j) 



(-1)* 



2 

i-iy 



4/^(2z)/^(2j) 

2?7 2^7 



Ji{2j) /i(2i) 

2?7 2?7 



(13) 



(14) 



where sgn(x) is the sign function and 
/a {x) = 



with 



2(L + 1) . / ir\x\ 
sin 



■ cot 



2(L + 1) V 2 ( L + 1) 



V2(L + 1) 

7T;r 



"2" = At) "T" — ^1! "2" = ^T-Z 



(15) 
(16) 

(17) 



(here and in the following all the distances are in units 
of the lattice constant). 

The agreement between exact numerical calculations of 
the XXZ correlation functions and analytical expressions 
in (13)-(14) is very good, and it becomes excellent with 
L ~ 100 for -0.8 < A < 0.8 43 . Thus one may readily 
assume that Eqs. (13)-(14) provide quite an accurate an- 
alytical expression for the spin-spin correlation functions 
in the XXZ model 84 . As a consequence, constructing a 
rigorous mapping between the BH and the XXZ spin-1/2 
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Hamiltonian and expressing correlation functions in one 
model in terms of the ones of the other gives an efficient 
and straightforward way to provide accurate analytic ex- 
pressions for real space correlation functions in the BH 
model at half-integer filling. 

We finally observe that the only system-dependent pa- 
rameter determining the spin-spin correlation functions 
is the coefficient rj: thus, in tracing out the mapping be- 
tween the two models, this is the key quantity to be calcu- 
lated as a function of the BH parameters. In particular, 
one may distinguish between the regions in parameter 
space with r) > 1/2 and r\ < 1/2: while the former one 
corresponds to an antifcrromagnctic spin chain, the latter 
one (which may be realized for pertinently chosen values 
of the parameters of Hbh , as we shall show below) cor- 
responds to a ferromagnetic chain. 



III. EFFECTIVE SPIN-1/2 HAMILTONIAN FOR 
THE BOSE-HUBBARD MODEL AT 
HALF-INTEGER FILLING 

As reviewed in the previous Section, for U — > oo, the 
BH Hamiltonian maps onto the XXZ model in Eq. (4), 
with the parameters J, A given in Eqs. (5, 6). This may 
be seen as a first-order term in an expansion (in powers 
of t/U) aimed at computing the effective Hamiltonian: 
in this Section we compute this effective Hamiltonian to 
the next order. As we shall show in the following, this is 
enough to fit quite well the numerical data for the cor- 
relation functions of the BH model using the analytical 
results obtained for the correlators of the XXZ chain. 

To approach the largc-JJ limit one may either proceed 
by performing a strong coupling expansion to the second- 
, or higher-order of perturbation theory, or by deriv- 
ing effective Hamiltonians using alternative techniques, 
based on canonical transformations or continuous uni- 
tary transformations 85 . Strong coupling techniques were 
extensively applied to the BH model at integer filling: 
as a result, one may, for instance, evaluate the energy of 
the Mott insulator and of the supcrfluid state in higher- 
order perturbation theory and determine the phase tran- 
sition point and the phase diagram in the U — fi plane 86 . 
Since we arc rather interested to the BH at half-integer 
filling, i.e., in the region of the phase diagram where 
the lobes touch and the superfluid phase persists also 
at very small U (with V = 0), we found it convenient 
to use an approach based on continuous unitary trans- 
formations 47 ' 48 . We follow the notation and the method 
presented in the paper by Glazek and Wilson (GW) 47 : 
systematically using the GW renormalization procedure, 
we work out an effective description of the dynamics of 
the BH model, restricted to the low-energy subspace de- 
termined by the constraint on the total number of par- 
ticles and by the large-?/ assumption. As a result, the 
low-energy subspace is spanned by states with cither n or 
n+1 particles per site, with the total number of particles 
being fixed to N. Thus, the space of physically relevant 



states at each site is in one-to-one correspondence with 
the Hilbert space of states of a quantum spin- 1/2 degree 
of freedom; we shall see that, at half-integer filling, even 
for finite U the BH model may be replaced by an effec- 
tive spin-1/2 Hamiltonian, with pertinently determined 
parameters. The method amounts to an iterative block- 
diagonalization of the BH Hamiltonian on the space of 
cigenfunctions of -Hbh with t = 0. 

To illustrate the procedure, we start from the explicit 
construction of the "low-energy" Hilbert space of physi- 
cally relevant states, in the large- U limit. Neglecting ex- 
citations with energy ~ U amounts to truncating the 
Hilbert space to a subspace J-", defined as 

T = Span{|n + fj,i, . . . ,n + hl)}, (18) 

with ^ taking the values fa =0,1 and Ylt=i Mi = 
In Eq. (18) \ni, ■ ■ ■ labels the state in the Hilbert 

space with m particles on site i. To implement the GW 
approach, one splits the Hamiltonian (3) as Hbh = -Ho + 
Hi, with 

H = -J2ni(ni-l) + Vj2 n i n i+i ( 19 ) 

i i 

Hi = -tY,(b\b i+1 + b\ +1 b t ). (20) 

i 

From Eqs. (19)-(20) one sees that H is diagonal with 
respect to the partition of the Hilbert space into J- plus 
its orthogonal complement, since 

H \m, ...,n L ) = E [m, . . .,n L ]\m,. ..,n L ) (21) 

with.E [ni,...,n L ] = (U/2) £\ rij(n,— 
while Hi exhibits off-diagonal (with respect to the par- 
tition of the Hilbert space) matrix elements which are 
0(tn). In order to block-diagonalize -Hbh, one needs to 
perform a similarity transformation 

Hbk -> Hbh = S^ BIJ S, (22) 

with S unitary. Upon setting S = I + T, the unitarity of 
S implies the optical theorem 

T + T t +T t T = 0. (23) 

Setting T = h + a, with 

h=i(T + Tt), a=i(T-Tt) (24) 

one finds that Eq. (23) yields 

h=l(a 2 -h 2 ). (25) 

Eq. (25) shows that h is always "higher order" than a. 
Following Ref. 47, it is most convenient to define the new 
interaction Hamiltonian Hi as 

Ui = H BH - H (26) 
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so that the new "free" Hamiltonian is the same as the 
old one (H ). 

To further proceed and determine S, one has to require 
that the matrix elements of Hi between states with en- 
ergy difference > U are equal to zero, amounting to state 
that T~Li is block-diagonal with respect to the partition of 
the Hilbcrt space into J- plus its orthogonal complement, 
i.e. 



VUiV + (l-V)Ui {l-V) = Hi, 



(27) 



where V is the projector onto T and I — V the projec- 
tor onto its complementary subspace 87 . One sees that 
Eq. (27) implies that 



THi {I-T) = {I- V) HiV = 0. 
Using Eqs. (22)-(24)-(26), one may write Hi as 



(28) 



Hi = (I + h - a) (H Q + Hi) (I + h + a) - H (29) 



and Eq. (28) then becomes 

V{Hi + {H Q , h} + [H , a] + T^Hi + H T T 
+ T t F / T} (I-V) =0. 

Eq. (30), together with the identity 

a = Pa(I-P) + (I-P)aP 



(30) 



(31) 



and with Eq. (25), is all what one needs in principle to 
fully determine a and h (and, therefore, the operator T). 

However, except for some simple cases 48 , an explicit 
solution for T cannot be exhibited. For this reason we 
proceed by writing the solution for T iteratively, in a 
series in Hi: in particular, we use Eq. (30) to determine 
a to first order (ai) in Hi. We provide the details in 
Appendix A and the result for ai in Eq. (A4). Using 
Eq. (A4) and setting T sa ai, we find that Eq. (22) reads 



S^bhS = F BH + [H ,ai] + [Hi, a x ]. 



(32) 



The GW procedure may be readily iterated to determine, 
in principle, T to any desired order in Hi. However, since 
keeping only second-order contributions in Hi provides 
already quite an excellent estimate for the real-space cor- 
relation functions of operators in the BH model (as ex- 
plicitly shown by the numerical calculations wc report 
in Section VI), setting T « a x already provides quite a 
good approximation to the exact T. 

Since the approach we are implementing is pertur- 
bative in Hi, one should enforce Eq. (30), as well 
as Eq. (31), to each order in Hi; moreover, since 
T'[Ho,ai]'P = 0, one may neglect the term [Ho,a.i] in 
Eq. (32) and approximate the effective Hamiltonian act- 
ing within T as 



The first term in the right hand side of Eq. (33) yields 
a spin-1/2 Hamiltonian which is actually the spin-1/2 
XXZ chain introduced in Section II B and having the 
anisotropy and the coupling given by Eqs. (5)-(6): 



H^ z = VH BU V = -JJ2^f +1 + s\s\ +1 - jsfsU) 

(34) 

with J = 2t (/ + ^) (constant terms have been omitted). 
The effective spin-1/2 operators are defined as 



v(-bi + bt)v, 



V[b\h-f)V; 



(35) 



the boson number eigenstates at site i correspond to 
the eigenstates of s\ according to \n)i <-> | and 
\n + l)i o | t)i- Therefore, the result in Eq. (34) cor- 
responds to the "naive" large-J7 limit for the BH model 
at half- integer filling discussed in Section II B, in which 
off-diagonal matrix elements of relevant operators (in- 
cluding the Hamiltonian itself) are set to zero from the 
very beginning. 

Corrections to -ffxxz arising from virtual transitions 
involving states outside of J- may be properly accounted 
for within GW procedure, allowing to get the effective 
spin-1/2 Hamiltonian to the next order in t/U. Summing 
over all virtual transitions outside of J- induced by Hi, 
one finds 



F« = T[Hi, ai ]T = -t 2 Y^V (b)b J+1 + b] +1 b 

3,1 



H cS = V {i/ BH + [Hi, ai]} V = H 



(o) 

XXZ 



H 



(i) 



(33) 



X (I - V) (Hbk)- 1 (I-V)[ b)b l+1 + b] +1 b t V. (36) 



In particular, when computing H^\ one has to consider 
intermediate states with either one of the fij in Eq. (18) 
being equal to 2, or to —1 (all these states have en- 
ergy ~ U, with respect to states in the subspace F), or 
states with one of the /ij equal to 2 (—1), and the other 
equal to —1 (2) (all these states have energy ~ 2U, with 
respect to states in the subspace F). Thus, one even- 
tually finds out that H^ can be written as the sum 
of two terms: H^ = H^ ag + H^ d , with H^ ag being 
the part of H^ having 1- and 2-nearest-neighbour spin 
terms, while H^g d contains 2-next-nearest-neighbour and 
3-spin terms. Omitting constant terms, their expression 
are given by: 
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H 



(i) 



H 



(i) 

offd 



4(n+ l)t 2 
U 

t 2 (n + l) 2 
U 



J2 Si - jj (3n 2 + 6n + 4) ]T 



XI ( s i+i s i-i + s *+i s ti) 



2t 2 (n+ 1) 

c7 



(37) 
(38) 



As we shall in the next Section, using a Luttingcr liquid 
representation, may be recasted in the XXZ form 

with coupling and anisotropy coefficients depending on 
U. 



IV. EFFECTIVE XXZ PARAMETERS VIA A 
LUTTINGER LIQUID REPRESENTATION 

The effective spin Hamiltonian in Eq. (33) is not in 
the XXZ form: in this Section we show how the con- 
tribution coming from H^ 1 ' may be accounted for by a 
pertinent redefinition of the parameters of the spin- 1/2 



XXZ-Hamiltonian H 



(0) 

xxz- 



The first contribution to -H^iag in the right-hand side of 
Eq. (37) describes an effective magnetic field in the z di- 
rection 88 , while the second term simply shifts the value 
of the XXZ anisotropy. At variance, the term H^ d in 
Eq. (38) contains 3-spin, as well as non-nearest neigh- 
bour, couplings. To show how these terms can be ac- 
counted for via a redefinition of -ffx°xz' ^ ^ s mos t con- 
venient to rewrite H^ d in terms of the Jordan- Wigncr 
(JW) fermions aj,aj, as outlined in Appendix B. As a 

result, one writes H^ d as a sum of a bilinear (H 2 ), plus 
a quartic (Hi) term, that is 



H 4 



with 
H 2 

Ha 



t 2 (n + l) 
U 

2t 2 (n+l)'- 
U 



X ( a Li«m + a t+i a *-i) ( 39 ) 

i 

(aj-i a i+i + al+i a *-i)( 4 °) 



Since H 2 is bilinear in the JW fermions, it merely mod- 
ifies the single-fcrmion dispersion relation, yielding the 
quadratic Hamiltonian in the JW fermions reading 



H 



(0) 
XXZ 



H 2 



E 



2 J cos k 



t 2 (n + l) 
U 



B }ala k 



cos (2k) 
(41) 

with B = A(n + l)t 2 /U. Setting e(k) = -2 J cos k + 
* cos (2k) — B, one finds that the Fermi points, de- 



fined by e(kp) = 0, are given by 



cos ki 



U(n + 1) fU(n + l) 



2,7 



V 2J 



+ n + 2. (42) 



Upon linearizing the dispersion relation around ±kp and 
setting k = kp + p, one gets 



e(±kp + p) ~ ±Jsinfci? 



1 - 



2 J 



U(n + 1) 



COS kp 



p. (43) 



From Eq. (43) one sees that, since cos kp 7^ 0, H2 implies 
a nonzero effective magnetic field B c s 8& , as well as a re- 
definition of the Fermi velocity vp. This yields a redefined 
coupling given by S cff /J cff = — cos kp. Since 



B c ff 
one obtains 



-J cos kp 1 



2J 
~U~ 



cos ki 



J, 



off 



2J 

J ( 1 — — cos kp 



(44) 



The quartic term H4 can be dealt with by noticing that, 
in the low-energy, long-wavelength limit, one can write 



a]_ x a j+1 + aj+iflj-i — > -^p R (xj) + p L (xj) 

-i-iYbi>U*i)i>L(xj) + , (45) 



where the chiral fermion fields ip R (xj), ifjL(xj) are de- 
fined from the long-wavelength expansion of dj as (see 
Appendix B) 



(46) 



with Xj = aj, and the chiral fermion densities given by 
Pr(xj) = ^ R (xj)ip R (xj) and pp(xj) = ijj' L (Xj)ip L (xj). As 
a result, Hi may be written as 



H 4 = - 



U 2 (n+lf 
U 



dx 



{(PB.(X)) 2 



+ (p L (x)) 2 +Ap R (x)p L (x)Y (47) 

Comparing Eq. (47) to Eq. (B13), one sees that H4 takes 
the same form as the term JY]j SjSj +1 in the spin-1/2 
XXZ Hamiltonian in Eq. (4). 
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1.5 - 



< 0.5 



-0.5 - 








n = 


— n = 


10 


n = 


+oo 




0.1 



0.2 

J/U 



0.3 



0.4 



From Fig. 1 one also sees that A e g- may be tuned by 
varying the ratio J/U: in particular A c ff can be differ- 
ent from even if V = (as it is typical for alkali 
atoms). Fig. 1 also suggests the possibility of describing 
the whole phase diagram of the XXZ spin-1/2 chain using 
BH model for a single species of bosons with pertinently 
chosen parameters, see also Section VI 89 . 

Finally we notice that, since the sign of A e g- may be 
changed by a pertinent choice of J/U and V, the Lut- 
tinger liquid effectively describing the XXZ-Hamiltonian 
may be repulsive or attractive. As noticed in the context 
of ID Josephson junction arrays 72,73 , the transition be- 
tween the repulsive and the attractive side may be moni- 
tored by inserting a weak link (i.e., a nonmagnetic impu- 
rity 82 ): it would be then interesting to analyze the effects 
of a weak link introduced in a bosonic system described 
by the BH Hamiltonian. 



FIG. 1: A ft vs. J/U for different values of V/J and n. The 
top (bottom) dotted line corresponds to the value of A e g for 
U/J -> oo and V/J = 0.5 (V/J = 0). The other lines are for 
V/J = 0.5 (top) and V/J = (bottom), with n = (solid 
black lines), n = 10 (dashed red lines) and n — > oo (dot- 
dashed blue blue lines). Inset: same as in the main panel, but 
for rj vs. J/U. 



Collecting together all the above results allows to write 
an effective XXZ Hamiltonian, describing the BH model 
to the order (t/U) 2 , as: 



rreff 
n XXZ 



V. CORRELATION FUNCTIONS 

The mapping between Hbh and Hxxz derived in Sec- 
tion IV enables to select the ground-states on which to 
compute the pertinent vacuum expectation values. In- 
deed if |$o) is the ground-state of the BH Hamiltonian 
given in Eq. (3), and |*o) = St|$o) is the ground-state 
of H e g = S^iJgjjS, the GW approach requires 

($ O |0BH [{M f }] $0> = (folS^BH [{M f }] S|* ) 

= (*o|Oxxz[K}]|*o), (52) 



with J e ff defined in Eq. (44) and 



ioff 



jf cos(fcp) 



with 



A 



V 
1 



t 2 (3n 2 + 6n + 4) 4t 2 (n + lf 



JU 



JU 



(48) 



(49) 



(50) 



Since J c s acts just as an effective over-all scale of 
Hxxz> then A e s is the only parameter determining the 
behavior of spin-spin correlations in the XXZ model. Sub- 
stituting Eq. (49) in Eq. (8) one gets 



1 

rj = 1 arccos A c ff, 

7T 



(51) 



which provides an explicit formula for the effective Lut- 
tinger parameter for the BH model at half-integer filling. 
In Fig. 1 we plot both A e s and r/ versus J/U, for dif- 
ferent values of V/J and n. One sees that fj, = 10 and 
?i — > oo are almost indistinguishable, and that the limit 
of the quantum phase model for Josephson junction ar- 
rays (n ^> 1) at offset charge q = e is practically reached 
at n ~ 10. Furthermore, one sees that the dependence of 
rj upon n is rather small. 



where Obh \{b, 6^}] (Oxxz [{s a }]) denotes a generic BH 
(XXZ) operator. Of course, Eq. (52) is exact only if S is 
the exact solution of the GW equation (31): by comput- 
ing it perturbativcly at a given order, one recovers the 
correspondence between ground-state expectation values 
of BH and spin-1/2 operators at the chosen order. 

In the rest of the paper, we will be interested in corre- 
lation functions of the following BH operators: 



M 



i -j 



(m - f) (nj - f) 
- hU 



(53) 
(54) 



Using the results of Appendix C one has S^A^f^S = 



M 



so that 



(*o| (ni - /) (n 3 - /) |* ) = (*o|*fa||*o) + O 



More generally, 
fies {1-V)O bb T> 
($o|0bh|$o) 



if the operator 
= Wbh (I - V) 



U 2 
(55) 

Obh satis- 
= 0, then 



(folCxxzl^o), with Oxxz obtained 
from Obh by substituting b\ and n, — / respectively 
with s7, sf and sf. At variance, for M.j~, one obtains a 
more involved expression (see Appendix C for details): 
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■oi&l&^o) « (n + i)(* |*r4i*o) + f(n + S ( " + 1) (^oik- + i + + [4+i + 4-i\ s ri*o> 



tn(n + 1) 
2£7 



I'/ 



2 + S5 



l*o>, 



where again we neglected contributions arising to 



VI. RESULTS 

In this Section we compare the numerical results ob- 
tained through DMRG for the correlation functions and 
the phase diagram of the BH model with the analytical 
predictions for the correlators from the effective Hamil- 
tonian i?xxz gi ven by Eq. (48). 



A. Correlation functions 

Let us focus on the BH correlation functions. Since 
DMRG simulations are made on a finite number of sites L 
and for open boundary conditions, we may use Eqs. (13) 
and (14) yielding the zz and xy correlation functions 
of the XXZ model. We evaluate the values of the non- 
universal constants a, b, c defined in Eq. (17) both nu- 
merically, by DMRG simulations of the XXZ model, and 
using the analytical expressions presented in Refs. 44,45 
and reported in Section II B. As confirmed in Ref. 43, the 
values of a, b, c obtained in the two ways are in excellent 
agreement. We show that the analytical expressions for 
the XXZ correlations are well confirmed by the numer- 
ical BH correlations also for small L (e.g., for L = 30) 
and for J/U relatively large (as large as ~ 0.5). It should 
be stressed that, at variance, the agreement is not very 
good by setting A c g = V/J, i.e. by using the Hamilto- 

nian H^xz obtained for U — > oo neglecting contributions 
arising from the GW procedure. 

The correlators ($o| {n% — f) (rij — /) |$o) and 
(&o\b\bj\$>o) are evaluated from the corresponding XXZ 
quantities using respectively Eqs. (55) and (56). They 
are plotted as a function of r = \i — j\, with i and j such 
that 43 i = (L — r + l)/2, j = (L + r + l)/2 for odd r, 
and i = (L — r)/2, j = (L + r) /2 for even r (for instance, 
for L = 100 sites, r = 1 corresponds to i = 50, j = 51; 
r = 2 corresponds to i = 49, j = 51; r = 3 corresponds 
to i = 49, j = 52, and so on). 

In Fig. 2 we plot our results for the density-density 
correlations ((rij — /) (rij — /)) for a typical set of val- 
ues, i.e. for U = lOt, V = 0.5i, / = 0.5, corresponding to 
J/U = 0.2. Black dots (joint by a black line as a guide 
for eye) are the density-density correlations evaluated in 
the BH model, red squares (line) are the correlation func- 



tions (sfSj) in the ground-state of the XXZ chain with 
effective anisotropy A e s given by (49) and the a, b, c con- 
stants numerically determined from DMRG simulations 
of the XXZ chain, while the green diamonds (line) cor- 
respond to a, b, c analytically determined from Eqs. (17) 
and (9)-(10)-(ll). One sees that, up to numerical accu- 
racy ^ 10~ 5 , the results obtained analytically for the 
XXZ effective model are in excellent agreement with re- 
sults of the density-density BH model even at small dis- 
tance. Blue triangles (line) display the XXZ Hamiltonian 
results in the U — > oo limit, with anisotropy A = V/J - 
one sees that the relative error is noticeably larger. 

In Fig. 3 we plot the off-diagonal correlations (b\bj) 
for the same set of values of the BH parameters (namely, 
U = lOt, V = 0M, f = 0.5, L = 150). The meaning of 
dots and lines is the same as in Fig. 2: also here one sees 
that the results obtained from the GW effective Hamilto- 
nian i?xxz are m m uch better agreement than the ones 
obtained using -f/xxz with A = V / J, this happens even 
though J/U is as low as 0.2. 

To quantify the agreement between BH and XXZ re- 
sults, we consider the absolute value of the relative error 
done in evaluating a correlator C(r) as the ground-state 
average of the corresponding operators in the BH model 
[Cbh(0], and in the XXZ model [C XX z(r)]. More pre- 




r 



FIG. 2: Density-density correlations \{(rii — f) (rij — /))| vs. 
r — \i ~ j\ f° r U — Wt, V = 0.5t and / = 0.5, with number of 
sites L = 150. Black dots and line: BH results; red squares and 
line: XXZ result with a, b, c numerically determined; green 
diamonds and line: XXZ result with a, b, c analytically de- 
termined; blue triangles and dotted line: U — s- oo XXZ result 
(indicated by the label "infmite-f7 limit" - see text fur further 
details). 
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150 



FIG. 3: Real part of {b\bj) vs. r = \i - j\ for {7 = Wt, V = 
0.5i, / = 0.5, L = 150. Magenta stars and line denote XXZ 
results with "non-rotated operators" . 




FIG. 4: Relative error of the correlation function 
\{(th — /) {rij — /)}| vs. r — \i — j\ for the same param- 
eters (and the same conventions for symbols and lines) of 
Fig. 1 - numerical XXZ red squares and analytical XXZ 
green diamonds for finite U almost coincide. The average 
value, with r max = 3L/5, is 0.06 ± 0.04 for the finite-?/ XXZ 
model and 2.9 ± 1.2 for the infinite-17 XXZ model. 



cisely, we define 



6C(r) = 



Cbh(V) - Cxxz(r) 



Cbh(^) 



(56) 



focusing on C zz (r) = \{(m - /) {rij - /))| and C xy (r) = 
Re &_,-)] . To summarize the information on the relative 
error, we compute the average value S^C and the stan- 
dard deviation of the relative error (56) for a distance 
r = \i — j\ between a minimum value r m j n = 1 (2) for zz 
(xy) correlations, and a maximum value r max ~ 3L/5. 

The relative errors for the zz and xy correlation func- 
tions are plotted in Figs. 4-5: the error made using the 
GW -ffxxz i s °f the order of few percents (in agreement 
with (J/U) 2 — 0.04). At variance, the relative error made 




FIG. 5: Relative error of the real part of {b\bj) vs. r — \i — j\ 
for the same parameters (and the same conventions for sym- 
bols and lines) of Fig. 2. The average values, with r max = 
3L/5, are: 0.025 ±0.004 (green diamonds - analytical XXZ re- 
sults), 0.38±0.08 (blue triangles - XXZ result in the infinite-C/ 
limit), 0.26 ± 0.03 (magenta stars - finite-17 result with non- 
rotated operators). 



by using the XXZ model in the infinite-[/ limit with- 
out applying the GW procedure is much larger, although 
the value of J/U is not so large. Indeed, the error (5 av is 
~ 300% for the zz correlations and ~ 40% for the xy cor- 
relations (to be compared with ~ 6% and ~ 3% obtained 
from i?xxz)- We checked that these results do not depend 
on the particular choice of r max : of course, when r max is 
closer to L, the error is larger (especially for the density- 
density correlations) due to boundary effects. From the 
data of Figs. 4-5, one also sees that, at short distance, 
it is larger than that at intermediate distances (with r 
being few units it is < 10%). As expected, it decreases 
at the center of the chain r ~ L/2, while, close to the 
end of the chain r ~ L, it increases. We also observe that 
finite-size effects are less visible for xy correlations. 

The agreement between numerical and analytical re- 
sults turns out to be stable also if one takes chains with 
smaller sizes, as it is apparent from Figs. 6-7, where we 
plot the zz and xy correlation functions for different sizes 
L. The corresponding errors are given in the following ta- 
ble: 



L 


8 {U) C 


u &v ^ZZ 


S (U) C 


x(«0 c 

u av i-'xy 


30 


0.12 ±0.08 


1.7 ±0.9 


0.04 ±0.01 


0.26 ±0.07 


50 


0.10 ±0.06 


2.0 ± 1.0 


0.04 ±0.01 


0.29 ±0.08 


80 


0.09 ±0.05 


2.4 ± 1.0 


0.04 ±0.01 


0.33 ±0.08 


100 


0.08 ±0.05 


2.5 ± 1.1 


0.04 ±0.01 


0.35 ±0.08 


150 


0.07 ±0.04 


2.9 ± 1.2 


0.04 ±0.01 


0.38 ±0.08 



where for simplicity sffi (#iv ) denotes the average er- 
ror for the XXZ correlators at finite-?/ (infinite-?/ limit) 
with (without) the GW procedure. We see that, for 
the density-density zz correlations, the average error in- 
creases when the size L decreases. 
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FIG. 6: In each panel we plot \{(rii — f) (rij — /))] vs. r — 
\i-j\ for different sizes: L = 30,50,80,100 (U = lOi, V = 
0.5t, f = 0.5). The blue data (filled triangles) denote the 
innnite-t/ XXZ results. 




FIG. 7: Real part of {b\bj) vs. r = \i — j\ for the sizes and the 
parameters of Fig. 6. 



In Figs. 3 and 5 wc also plotted (magenta stars and 
lines) the results obtained according to Eq. (52), where 
we took Oxxz = Cbh and not Oxxz = S^ObhS. Indeed, 
as stressed in Ref. 90, solving the equation for S amounts 
to perturbatively find a transformation enabling to block- 
diagonalize i?BH- The ground-state of Hbh changes ac- 
cordingly: if one wants to compute expectation values 
of certain operators in the BH model, one has to rotate 
the chosen operator according the S transformation - in 
other words, physical quantities in the effective theory 
are not simply the expectation values of the operators in 
the projected subspace: this guarantees the unitarity of 
the procedure. An example is already provided in Ref. 50 
for the computation of the staggered magnetization in 
the 2D Fcrmi-Hubbard model with large- U effective spin 
models. In Appendix C we give details on the the explicit 
computation of the GW rotation for the operators b[bj 
and (n, — /) (rij — /). We remark that, while for density- 
density correlation functions (zz correlations in the XXZ 
model) magenta stars coincide with the black circles, this 
is not the case for (b\bj) (xy planar correlations in the 



FIG. 8: Density-density correlations \{(rii — f) (rij — /))| vs. 
r = \i — j\ for different values of U/t = 20,10,5,3.3, cor- 
responding, respectively, to J/U = 0.1,0.2,0.4,0.6 (with 
V = 0.5t, / = 0.5, L = 150). The infinite-!/ result in each 
plot is the top blue line (triangles), while the other two - al- 
most coinciding - lines are the numerical BH (black circles) 
and XXZ (red squares) results. 



U=20t 



U= 10 t 




so r ioo 



so r ioo 



FIG. 9: Real part of {b\bj) vs. r = \i — j\ for the same pa- 
rameters (and the same conventions for symbols and lines) of 
Fig. 8. 



XXZ model). 

In Figs. 8-9 we plot the zz and xy correlation func- 
tions for different values of U: in these plots the ratio 
J/U ranges from 0.1 to 0.6. As expected, one sees that 
for J/U = 0.1 the relative error made by the infinite-C/ 
results is not very large (~ 10% for (b\bj) correlations), 
but, as soon as J / U > 0.2, it is already well visible. The 
relative error made by using the effective -ff X xz turns out 
rather small even for J/U = 0.6, where the relative error 
on density-density correlations is only sa 7%, while for 



(b\bj) correlations it is w 15% 
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B. Antiferromagnet and domain ferromagnet in 
the ID Bose-Hubbard model 

The XXZ model is gapless and critical for — 1 < A < 
1, antifcrromagnctic for A > 1 and ferromagnetic for 
A < — 1: in the latter ferromagnetic phase, all the spins 
are aligned. However, the BH at half-integer filling maps 
into the effective XXZ chain (48) supplemented by the 
condition that the total spin is vanishing: therefore we 
expect that, in the BH model at A e g- < — 1, domain walls 
form separating regions with "up" spins (i.e, with f + 1/2 
particles per site) and regions with "down" spins (i.e., 
with / — 1/2 particles per site). At variance, at A e ff > 
1 the staggered magnetization becomes non vanishing: 
in the bosonic BH language the antiferromagnetic state 
corresponds to the "charge-checkerboard ordered state" 
|/ + 1/2,/ -1/2,/ +1/2,/ -1/2,..-). 

This shows that, consistently with the XXZ represen- 
tation of the BH model at half-integer filling, at finite U 
one can realize the transition between the spin-liquid and 
the Neel-Ising antiferromagnetic phase of the XXZ model 
(superfluid-to-charge density wave phase transition of the 
BH model), as well as the transition between the spin- 
liquid and the domain ferromagnetic Ising phase of the 
XXZ model (superfluid-to-domain Mott-insulating phase 
transition of the BH model) 89 . Since the former transition 
sets in at A c g = 1 and the latter one at A e g = — 1, using 
Eq. (49) for A e ff allows to determine the corresponding 
phase boundaries in terms of the parameters of the BH 
Hamiltonian. 

A complete discussion of the phase diagram of the BH 
chain in presence of nearest-neighbour interactions is pro- 
vided in Ref. 63: here we just focus on the half-integer 
BH chain with parameters chosen so as to lie close to 
A c ff = ±1, in order to show that the effective XXZ rep- 
resentation given in this paper also provides a good de- 
scription of these transitions. 

For the spin liquid-ferromagnetic transition, we stud- 
ied the BH chain with open boundary conditions vary- 
ing V (similar results are obtained varying t) and we 
plot in Fig. 10 the expectation value of (n, — /) as a 
function of the position i along the chain. We observe 
that, as a consequence of the open boundary conditions, 
a magnetic field proportional to V on the two bound- 
aries (i.e., at i = 1 and i = L) appears, whose effect close 
to the boundaries is clearly visible in the figure. Com- 
puting the quantity M = E r (-l)'~ J ((n, - /)(n,- - /)), 
one sees that it significantly increases around a critical 
value A^jT. From the numerical data for the BH model 
shown in Fig. 10 one may estimate A^jT ~ 1.05, in good 
agreement with the analytical value A^T = l 95 . We no- 
tice that a better estimate of A^ could be performed 
by adding a magnetic field in the boundaries to compen- 
sate the boundary magnetic fields arising from the open 
boundary conditions. 

Regarding the domain ferromagnet-superfluid transi- 
tion, we performed numerical simulations on the BH 
model with parameters chosen such that A c ff is close to 




I 



FIG. 10: Plot of ( — l) , + 1 (rii — /) vs. i numerically computed in 
the BH chain for U = Wt, f = 0.5 and L = 150 for different 
values of V: from top to bottom V/t — 3.3 (black circles), 
V/t — 3.1 (red squares), V/t = 2.9 (green diamonds) and 
V/t = 2 (blue stars) corresponding, respectively, to A e g = 
1.09,1.00,0.91,0.52. 

— 1 (sec Figs. 11-12). In Fig. 11 we plot (rii — f) as a 
function of the position i: one sees that the expectation 
value of the spin is constant and it changes sign close 
to the edges of the chain in order to satisfy the con- 
straint on the number conservation. For this reason we 
then plot the modulus of the same quantity in Fig. 12: 
since the average of the sf expectation values is of course 
zero, to determine the transition point from BH numeri- 
cal data we consider the averaged quantity X)j=i I ( n i~.f) I 
(e.g., for the different values of V shown in Fig. 12, such 
quantity is reported in the caption). From these data one 
can estimate that the domain ferromagnet is occurring at 
A^j ~ —1.02, with an error of few percent with respect 
to the analytical result A^ = — l 95 . Notice that the er- 
ror made by using -ffxxz with A = V/ J in the infinite-C/ 
limit is w 20%: as expected, the errors made in using the 
infinite- U results are generally smaller when one deals 
with global quantities. 



VII. CONCLUDING REMARKS 

In this paper we studied an XXZ representation of the 
Bose-Hubbard chain at half-integer filling for finite on- 
site interaction energy U . The effective XXZ model is 
obtained in two steps: first, we used a similarity renor- 
malization group procedure amounting to solve perturba- 
tively up to the order (t/U) 2 the exact equation for the 
operator block-diagonalizing the Bose-Hubbard model. 
The resulting spin-1/2 effective Hamiltonian is then re- 
casted as a XXZ spin-1/2 Hamiltonian with pertinently 
redefined coupling and anisotropy parameters. 

We use this mapping to provide analytical estimates 
of the correlation functions of the Bose-Hubbard model 
at half-integer filling and finite U. We then compared 



14 




FIG. 11: Plot of (rii — f) vs. i numerically computed in the BH 
chain for U = Wt, f = 0.5 and L = 150 for different values of 
V: V/t = —1.3, —1.5, —1.6, —1.7, corresponding, respectively, 
to Aotf = -0.91, -1.00, -1.04, -1.08. 



V=-1.3t 



V=-1.6t 




FIG. 12: Plot of the modulus of (rii — f) vs. i computed in 
the BH chain for the same parameters of Fig. 11: V/t = 
— 1.3, —1.5, —1.6, —1.7. The corresponding averages are w 
0.02,0.04,0.06,0.21. 



these analytical results with the outcomes of the numer- 
ical DMRG evaluation of the Bose-Hubbard correlation 
functions. We found that the agreement is very good, 
also for J/U rather large (~ 0.5) and for small num- 
ber of sizes (L ~ 30). Such a good agreement is not 
achieved - even for J/U relatively small (~ 0.1) - if one 

uses the XXZ Hamiltonian ff£xz with J = 2t(f + 1/2) 
and A = V/ J corresponding to the infinite-coupling limit 
of the Bose-Hubbard model. The transitions predicted at 
A e fi = ±1 for the XXZ chain are as well compared with 
Bose-Hubbard results, and a good agreement is found. 

Since the BH model at half-integer filing is not inte- 
grable or exactly solvable, it is quite valuable to have 
analytical estimates for its correlation functions. Besides 
its mathematical interest, we stress out that our results 
can be viewed from a two-fold point of view: on one side, 
we use known results from the (integrable) XXZ model 
to construct with high accuracy correlation functions of 



the Bose-Hubbard model. On the other side, the Bose- 
Hubbard chain at half-filling and at finite U may be seen 
as a quantum simulator of the XXZ chain. 

In our approach, the effect of an harmonic trap re- 
sults in a locally varying magnetic field: we feel that it 
would be interesting to compare the results stemming 
from an XXZ-based approach with the ones known in lit- 
erature for hard- and soft-core bosons in harmonic traps 
in the scaling limit 92 . In this paper we focused on the 
half-integer filling Bose-Hubbard model, but deviations 
from such filling could be easily accounted with the in- 
troduction of a magnetic field. We stress that the similar- 
ity Hamiltonian renormalization procedure could also be 
applied to bosonic ladders 93 and at integer filling, where 
a spin-1 model is found in the infinite- U limit. 

The large- V effects of edge magnetic field could also be 
studied, following the results known for the XXZ chain 94 : 
we observe that, for open boundary conditions and finite 
V, two boundary magnetic field terms — i3(,(sf +s|J, with 
Bt, oc V, emerge in the XXZ effective Hamiltonian 96 . 
Since a magnetic field at the edge induces corrections 
to the average value of sf decreasing as a power law 94 , 
these corrections are not only expected, but could be also 
worth the effort of future investigation. 
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Appendix A: Perturbative solution of the GW 
equation 

In this Appendix we show how use Eq. (30) to deter- 
mine a to first order in Hi, that is ai. To this order, one 
gets 



V{H I + [H ,a 1 }}(I-V) = 0, 
which may be solved by setting 



(Al) 



Pai (I - V) = 

VH! (I - V) { - TH T + (I - V) H (I - V) } _1 
+ [PH V,V ai (I-V)] 

k{-VH V+(I-V)H (I-V)}~ 1 . (A2) 

Up to term that are second order in tn/U, we 
may make the approximation VHqV « £o[^]Ij 
with £o[fh] = L {^n(n — 1) + Vfi 2 }, which implies 
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[PH V, Ta (I - V)] = 0. As a result, we get 
V&x (I - V) = VHj (I - V) 

x{-vH r + (i-r)H (i-r)}- 1 . (A3) 

Using the fact that a is antihermitean, from Eq. (A3) one 
obtains 

a x = VHi (I - V) 

x{ -VH V + {l-V)H {l-V)Y 1 - 

{-VH V + (I-V)H (I-V)y 1 
x(I-V)H I V. (A4) 



Appendix B: Jordan- Wigner fermion representation 
for the spin-1/2 XXZ spin chain 

In this Appendix we review the basic ingredients for 
fcrmionizing the spin-1/2 XXZ chain by means of JW 
fcrmions, which is used in the main text derive the effec- 
tive parameters to i?xxz ( see Section IV). 

We start from the spin-1/2 XXZ Hamiltonian with pe- 
riodic boundary conditions: 

L 

#xxz = ~ J { s t s i+i + 4+i s i - As " s m) - Bs l 

i=l 

(Bl) 

where s"_i_x = s°. The JW-fermion representation is re- 
alized in terms of lattice complex fermions ai,a\, 

s + = a y«i:iz\aUr } sf = a t a ._i. (B2) 

The first step is to fcrmionizc the s + s~ term in the stan- 
dard XXZ spin-1/2 Hamiltonian plus a possible magnetic 
field contribution: this yields ifxxz — > H% with 

H 2 ] = E - J ( a l a ^ + 4+1°*) - B - I) ■ 



Defining the Fourier modes as 

L 

J k 3, 



= 7l ^ e " 



3=1 



one gets 



(B3) 
(B4) 

(B5) 



with eo(fc) = —2Jcosk — B. Since the effective XXZ spin- 
1/2 Hamiltonian we derive in this paper is the effective 
description of a BH Hamiltonian at fixed filling in the 
canonical ensemble, the z-componcnt of the total spin is 
set to and the magnetic field term does not provide 
any contribution to the total Hamiltonian. Accordingly, 



we may set B = from the beginning. This implies that 
the Fermi momenta are given by ±/cf = where a is 
the lattice constant. As a consequence, the effective low- 
energy, long-wavelength Hamiltonian may be obtained 
by expanding the lattice fermion fields around the Fermi 
points ±fci? = ±5, as 



{e ikF ^ R ( Xj ) + e- ik ^Mxj)}, 



(B6) 



with 



^L{Xj) 



bl<A 



-y 



sfL 



(B7) 



bl<A 



Xj being equal to aj, while A is a pertinent cut- 
off in momentum space. Accordingly, the low-energy, 
long-wavclcngth quadratic part of the fermionized XXZ 
Hamiltonian will be given by 



H. 



(0) 



ivp I dx < if> L (x) 



dx 



dx 



(B8) 

with the Fermi velocity Vp = 2aJ. Moreover, since the 
ground-state |G) of ^ s defined by 

|G)= J] 4|0>, (B9) 

co(fe)<0 

one gets that 



t ^ t 

a j a 3 - 2 J =: a 3 a 3 : ' 



(BIO) 



where the colons denote normal ordering with respect to 
|G). Accordingly, one may rewrite the "Ising" contribu- 
tion to the spin-1/2 XXZ Hamiltonian as 



L-l L-l 

fflsing EE JA S i S l+l > JA Yl 'A* 1 !'- - a )+l a 3 + l- ■ 

i=l i=l 

(Bll) 

Resorting to continuum chiral fermion operators, from 
Eq. (B7) one obtains 



:a)a 3 : 



:^ R (x 3 )^ R (x 3 ): + :if)l(xj)i/)L(xj): 



Setting 



:a] aj : 



p R {xj) + PL {x 3 ) (B12) 
+ (-!) J {^U x 3)Mx 3 ) + ^{(x 3 )ib R (x 3 )} 
and summing over j one eventually gets 

L 

Rising ~> JAa / dx 3 {(p R (x 3 )) 2 + (pl(xj)) 2 



+Ap R (x j )p L {Xj)} + -Humklapp- (B13) 
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The term -Humkiapp in Eq. (B13) arises from two-fermion 
scattering processes in which the total momentum ex- 
changed by the particles is ~ 4fcp = — . A systematic 
analysis of such a term may be provided within the frame- 
work of a Luttinger liquid theory, obtained by bosoniza- 
tion of JW fermions 29 . As our work is mainly concerned 
with the correspondence between the BH and the XXZ 
Hamiltonian, here we just point out that the Umklapp 
term may be either irrelevant, or relevant, according to 
whether A< 1/2, or A > 1/2. As soon as it becomes rel- 
evant, it triggers the phase transition between the spin 
liquid and the antiferromagnetic Neel phase of the XXZ 
model which, as we discussed in the main text, corre- 
sponds to the transition between the superfluid and the 
charge density wave phases of the BH model. 



Appendix C: GW transformation of operators 



An advantage of the GW procedure is that it may be 
easily applied to single-boson operators: in particular, we 
are interested in the average values of the operators M.^ 
and M-i t j defined in Eqs. (53) and (54). Since ai is fully 
off-diagonal and VMf j (I — V) =0, if one approximates 
T with ai, one obtains S^A-lf jS = M.\y Instead, acting 
onto M-ij gives raise to a more complicated expression: 
expressing the final result in terms of spin-1/2 variables, 
one obtains 



J 



tn{n + 2) (1 



U 



t(n + 2)(n+l) 
U 



2 



i+i 



t(n + 2)(n + 1) f _ , fl 
+ ^ 



2 ' « 



»i-l a t 1 



1 



3+1 \ 2 

r 



i 

2 

< u 



(CI) 



2 J 



We observe that due to the constraint on the fixed to- 
tal particle number N, the total magnetization in any of 



the physical states of iJ e ff = -^xxz + ^d^ag + -^offd ^ s zcr0: 
since H e g contains no terms breaking the parity symme- 
try (sf —> —sf), its ground-state \^f ) is nondegenerate 
and, thus, it must be parity invariant. As a consequence, 



r(i) 



r (i) 



the average of any product of three spin-1/2 operators 
must necessarily give 0, greatly simplifying the calcula- 
tion of the ground-state average of the operator. 

Using this result one can obtain a simplified expression 
for ($o|A4^-|$o) at O onc has 



J 



(*o|^|*o>«<J|i_i|,i<*o| Q 
, i(n + 2)(n + l) 



U 

tn(n + 1) 



i+i M o + s * I 



S * S J + 1 { 2 + S J 



(C2) 



r 



Since any product of three spin-1/2 operators must nec- the main text, 
essarily give 0, then Eq. (C2) gives Eq. (56) reported in 
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